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We study the statistical and dynamical behavior of turbulent Kelvin waves propagating on quan- 
tized vortices in superfluids, and address the controversy concerning the energy spectrum that is 
associated with these excitations. Finding the correct energy spectrum is important because Kelvin 
waves play a major role in the dissipation of energy in superfluid turbulence at near-zero temper- 
atures. In this paper, we show analytically that the solution proposed in Ref. [l[ enjoys existence, 
uniqueness and regularity of the pre-factor. Furthermore, we present numerical results of the dy- 
namical equation that describes to leading order the non-local regime of the Kelvin wave dynamics. 
We compare our findings with the analytical results from the proposed local and non-local theories 
for Kelvin wave dynamics and show an agreement with the non-local predictions. Accordingly, the 
spectrum proposed in Ref. [l[ should be used in future theories of quantum turbulence. Finally, 
for weaker wave forcing we observe an intermittent behavior of the wave spectrum with a fluctuat- 
ing dissipative scale, which we interpreted as a flnite-size effect characteristic to mesoscopic wave 
turbulence. 

PACS numbers: 67.25.dk,47.37.-|-q,67.10.Fj 



I. INTRODUCTION 



When sufficiently large waves propagate through a 
medium, non-linear interactions cause an energy transfer 
between different scales; the energy that is introduced 
typically in a narrow band of wavelength is redistributed 
to waves of different wavelengths. Starting in the 1960s 
with oceanographers and meteorologists, 'weak wave tur- 
bulence' has now become a standard theory to explain 
the dynamical and statistical properties of an ensemble 
of weakly nonlinear interacting waves . Thanks to the 
generality of its formulation, this theory is an interdis- 
ciplinary subject; it has been implemented in situations 
ranging from Alfven waves in solar winds, ocean swells 
on a stormy sea, quantum waves in Bose-Einstein con- 
densates, Rossby waves in the atmospheres of rotating 
planets as well as in the thunder-like sound produced by 
vibrating thin elastic sheets, etc. Up-to-date discussions 
of the current state of wave turbulence are available for 
example in a recent book and in reviews [J, Q ; these 
sources also contain extensive lists of references. 

Quantum turbulence is the study of turbulent be- 
havior in zero temperature superfluids such as helium 
II d, [H, [l^ . One of the key properties of quantum tur- 
bulence is that it comprises of an inviscid flow which 
permits energy to reach scales far smaller than achieved 
in classical turbulence. At large scales, quantum turbu- 
lence shows a great similarity to classical turbulence - 
with the polarization of vortex bundles acting like large 
scale eddies. However, at small scales the analogy begins 
to break down. In the absence of viscosity in superflu- 
ids and at very low temperatures where thermal excita- 
tions (mutual friction) can be neglected, there is no clear 
dissipation mechanism that cuts off the energy cascade 
to very small lengthscales. Therefore, unlike in classical 



turbulence where viscous effects dissipate the energy, it 
is believed that for superfluids the energy is transferred 
to propagating waves on quantized vortex lines that are 
forced by vortex reconnections. Note that these so-called 
"Kelvin waves" were originally introduced in 1880 in the 
context of the classical Rankine vortex model [l3|- The 
Kelvin waves then interact allowing for energy to be 
transferred to higher frequency Kelvin waves until the 
wave frequency is sufficient for the excitation of phonons 
and thus the degradation of energy into heat [isj . 



Because of its importance in superfluid turbulence and 
the growing experimental capabilities in this field [M-Q, 
there has recently been a renewed interest in the statisti- 
cal physics of waves propagating on a vortex line 
A complete understanding of the statistical behavior of 
Kelvin waves is therefore crucial in order to develop a 
theory of superfluid turbulence. Until recently, Kelvin 
wave interactions were thought to be a local process - 
only Kelvin waves of similar wave numbers interact with 
each other thus forming the range of a local energy cas- 
cade 



11 



^ ^ The Kelvin wave cascade was described us- 
ing wave turbulence theory 0, Q for weakly interact- 
ing waves, resulting in the formation of a kinetic wave 
equation describing 3 O 3 Kelvin wave interactions [ll[. 
However in Ref. [12| . the locality assumption used in 
Ref. ll| was checked and was shown to be violated for 
local 3 3 Kelvin wave interactions. This invalidated 
the local theory, and a non-local theory was proposed [l| 
resulting in 1 -f-?- 3 Kelvin wave interactions. This has 
prompted a lively debate about the correct spectrum of 
Kelvin waves [m [13, [H [13, [13 . 
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A. Physical background 



The first term 



It is now recognized tliat the typical turbulent state of 
a superfluid consists of a complex tangle of quantized vor- 



equation 
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This equation contains a singularity as s — >■ r, and hence 
the model are postulated by introducing a cut-off in the 
integration: oq < |t" — s|, where on is interpreted as the 
core radius of the vortices [l^, [l^- Here, k — 2Trh/M 
is the quantum of circulation with M being , mass of 
''He atom (if the superfluid is '^He then M is the mass of 
two atoms, 27713). 

An important ste p in studying Kelvin wave turbulence 
was done by Sonin j2^ and later by Svistunov (23| who 
found a Hamiltonian form of the Biot-Savart Eq. ([T]) for 
a straight vortex line aligned in the z direction, e.g. line 
{x,y) = (00). Perturbing the line by small disturbances 
in the (x, ?/)-plane. 



w{z, t) — x{z, t) + iy{z, t), 



one writes: 
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where 6{. . . )/5w* is the functional derivative of (. . . ) and 
the superscript " * " denotes complex conjugation. The 
Hamiltonian for the Biot-Savart equation H is the energy 
of the system 
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Here we have used the notation w'(z) = dw/dz. 

Like in most of the examples of wave turbulence 
mentioned in the introduction, the theory for Kelvin 
waves starts by writing down a Hamiltonian equation for 
the complex canonical wave amplitudes af.(t) and a^(t) 
which are classical analogs of the creation and annihila- 
tion Bose operators in quantum mechanics [igI [iTj . 

Consider an isolated straight vortex line on a periodic 
domain of length C. One can write 0, Q 



dakjt) 
dt 
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where w{z,t) ~ K^^/-^^j^ak{t)exp{ik ■ z). The new 
Hamiltonian H{a,a*} = H{w,'w*} / C is the density of 
the old one and is a function of all ak{t) and a^(i) taken 
at the same time. 

For small Kelvin wave amplitudes (inclination angles) 
the Hamiltonian can be expanded with respect of , a^. 
The explicit form of iJ, Eq. (jSj, dictates an expansion of 
% in flfc and with even powers only, 



T~L — 'H,2 + T~L'\. 
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describes free propagation of Kelvin waves with a fre- 
quency Wfc. In turn, should be expanded in inverse 
powers of the large parameter A: 
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A = ln(£/a) 



(5c) 



Here a is the vortex line diameter and £ is the mean inter- 
vortex distance at which the description of Kelvin waves 
propagating along an individual vortex line fails. In typ- 
ical experiments A, in both '^He and "'He, is between 12 
and 15 @. It can be shown 11 1 that the leading approx- 
imation in A gives no energy exchange between Kelvin 
waves and therefore one has to account in Ji for sublead- 
ing terms, zero order in A, denoted by the superscript 

The higher order expansion terms in "Hint, ^4 and Hq, 
describe 2 <-> 2 and 3 o 3 scattering of Kelvin waves: 



J X! Tl't aia2alal , Uj = a{kj,t) , {5d) 



1+2=3+4 



^6 " ^ ^ifi aia2a3alalal , (5e) 

1+2+3=4+5+6 

Equations for the terms of order A^ and A'' in the inter- 
action amplitudes = r(fei, fc2|fe3, ^4) and W^i'2'3 = 
W{ki,k2,k3\k4,k5,k(i) (denoted as ^tI''^, ^tI^^ and 
^W^i'l, ^M^t'lJ.) were found in Ref. ^ and later con- 



firmed in Ref. PJ. 

On the face of it, the leading term in Hint (i-C ^4) 
describes a (2 ++ 2) scattering; the subleasing term, He, 
is responsible for the (3 ++ 3) scattering. It was argued, 
however, that the (2 ++ 2) scattering cannot redistribute 
energy between different scales. Therefore the (3 ++ 3) 
scattering becomes the leading interaction responsible for 
the inter-scale energy transfers. This does not mean, 
however, that the H4 Hamiltonian can be cornpletely 
disregarded. Instead, as explained in Refs. 0, [l3|, in- 
effective Hamiltonian (in our case Jii) can be eliminated 
from the problem by a proper non-linear canonical trans- 
formation {a, a*} ^ {b, b*}. This comes at a price of the 
appearance of "correction terms" in the 6-wave mixing 
amplitude (which we have denoted with a calligraphic 
vertex W, different from W) in the Hamiltonian He'- 

a = ^i^kbkbl+ne, (6a) 

k 

^6 = ^ E ^'1:2,3 bib2b,blblb; . (6b) 



1+2+3=4+5+6 



Note that the 4-wave mixing Hamiltonian has been 
eliminated altogether by the canonical transformation 
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{a, a*} ^ {b,b*}. The remnant of this Hamihonian ap- 
pears as an additional contribution to W. There is an 
exact relation between these two vertices which is best 
represented by a schematic graphic notation as 



Here we introduce the energy density in the fc-space, 



+ 72 



yy4,5,6 



72 { T^/n } . (6c) 



The 72 additional contributions arc schematically indi- 
cated as 72{...} since the exact expression is too long 
to be written here (one can handle them with a sym- 
bolic computation software such as Mathematica) . One 
observes that, on the one hand, they are 6-wave mixing 
terms, but due to their internal structure they can also be 
understood as the pairs of 4- wave mixing amplitudes (me- 
diated by a Green's function 1/fi). In the graph shown 
in (j6cp . with incoming wave vectors fci, ^2, starting 
from above and outgoing wave vectors k^, k^, fcg, the fre- 
quency = La{ki)+ijj{k2) — ijo{k4) — La{ki+k2 — fc4). The 
72 contributions of the type shown in ((6c)) differ by the 
directions of the arrows and by relabeling in fci, fc2, ^3 
and ^4, ^5, fee groups. 

This understanding provides us with an effective 6- 
wave mixing, allowing one to use standard procedure 0, 
0] for the statistical description of weak turbulence of 
Kelvin waves. It is based on the assumption that the 
turbulent dynamics of waves with small enough ampli- 
tude is chaotic and creates its own ergodic measure. The 
simplifying nature of weak wave turbulence is that, due 
to the existence of a small parameter, the statistical de- 
scription closes upon itself in terms of the pair correlation 
function 



(7a) 



which is also called the "wave action". Hereafter the 
pointed brackets stand for an average over the ergodic 
measure. The aim of the theory is to analyze the solu- 
tions of the kinetic equation of motion which is typically 
expressed as 



dukit) 
dt 



St{k,{nk'it)}). 
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where term St(fe, {rife/ (t)}) is so-called collision integral 
with an integrand proportional to the square of the ef- 
fective interaction amplitude W. This term is a function 
of k and a functional of rifc/ {t) 0, Q . 

The main part of the Kelvin wave energy E in the 
regime of weak wave turbulence is given by E2 = {'H2), 
where H2 is defined by Eq. ([5b)) . Together with Eq. (|7a)) 
this gives 
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E{k) = Ek^ LJkNk ., 
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which traditionally is called "energy spectrum" and we 
define the "wave action spectrum" as Nk = nk + ?i-fc, 
where k = \k\. According to Eq. ([5c)) in the leading in 
A approximation the energy spectrum of Kelvin waves is 
related to the wave action Nk as follows: 



Ek = ^k^Nk 



(8c) 



Up to this point these considerations are agreed by one 
and all, and are the basis of further developments. 



B. The controversy 

Recently Kosik and Svistunov derived an energy 
spectrum of Kelvin waves turbulence: 



E^,ik)=C^ 
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where e is the energy flux over scales and C^g is 
yet unknown dimensionless constant. Later L'vov and 
Nazarenko derived a very different result for the same 
spectrum: 



E^^ik)=C, 
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Here C^^ is another dimensionless constant. Both spec- 
tra are supposed to be "universal" , (i.e. independent of 
details of the energy forcing) in the the scaling range 
fc > fcf , where fc^ is the forcing wave number. 

The disagreement between the spectra ^ resulted in a 
heated debate concerning the correct nature of the energy 
spectrum. To identify the origin of this controversy we 
should clearly state that both results ([9]) were obtained 
within the same formal setup described in the previous 
subsection under the same set of assumptions about small 
nonlinearity and random phases. Therefore the difference 
between the spectra must originate from one or more 
mistakes made by either or both derivations. Indeed, we 
will argue that the mistake leading to the wrong result 
(|9al) is in the wrong assumption about the asymptotic 
behavior of the effective amplitude Wi'2'3 in the region 
where one of the wave vectors in the interacting sextet is 
much smaller than at least one other wave vector from 
the same sextet. The form of this asymptotics crucially 
affects the nature of the energy transfer in the Kelvin 
wave cascade. 

The derivation of the interaction vertex Wi'2'3 for 
Kelvin wave turbulence is not an easy task. In Ref. [ll| 
the explicit form of W^'2'3 was not presented. Instead, 
the authors "simulated the collision integral by a Monte 
Carlo method" with the conclusion that it converges and 
the main contribution to the energy evolution of Kelvin 
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wave with given fc-wave vector originates from the en- 
ergy exchange with other Kelvin wave with k' ^ k. This 
statement of the locality of the energy transfer allowed 
the authors to use a dimensional estimate that leads to 
the the spectrum ([9a|) . 

This result was criticized in [l2| where an explicit 
expression for the interaction amplitude yVi'2'3 (in the 
asymptotic region) was derived: 



3,4,5 

fc,l,2 
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Based on this equation, it was shown analytically in 
Ref. that the collision integral diverges. Moreover, as 
found in Ref. [l^ . two important sets of terms in 
of the order of unity were overlooked in Ref. [llj. One 
was the consequence of a trivial algebraic mistake in the 
Taylor expansion in Eq. (jlOp . where the authors forgot 
to expand the denominator. In addition, in Ref. 
there was a conceptual mistake: in the kinetic equation 
(see below) the conservation of energy requires the ex- 
act frequency-resonance condition that accounts for the 
sub- leading contribution ^ujk, ([5c|. while cancellation of 
linear in A terms in W^'^'j takes place only on the local 
induction approximation manifold. Remaining contribu- 
tions of the order of unity to where omitted in 
Ref. It remains unclear, how Ref. [ll| succeeded to 
state convergence of the collision integral; this must be 
either due to a mistakes in the calculation of Wfc'f'2 or 
due to an inaccurate implementation of the Monte Carlo 
numerical procedure. 

The divergence of the collision integral makes the as- 
sumption of locality for the KS spectrum invalid p^ . 
The non-locality of the spectrum (|9a|) implies that it is 
un-realizable and physically irrelevant. Recognizing the 
need for a greater understanding of the correct interac- 
tion term, a new non-local theory was promptly put for- 
ward, suggesting the alternative energy spectrum (j9b[) of 
Kelvin waves. 

This result was immediately attacked in Ref. psl ]. 
based on some nai've symmetry arguments claiming that 



(10b) 



for small k instead of the linear asymptotics (|10a|) . Note 
that in the case of quadratic asymptotics (|10b|) the colli- 
sion integral converges. Then the Kelvin wave energy cas- 
cade would be dominated by local interactions and spec- 
trum (|9ap would be valid. In the case of linear asymp- 
totics (|10ap the collision integral diverges [l^l , the Kelvin 
wave energy cascade is dominated by different-scale in- 
teractions and the spectrum ((9b| takes place [l[ . 

To conclude, the controversy regarding the Kelvin 
wave energy spectra ([9]) is a direct consequence of the dis- 
agreement regarding the asymptotic behavior of W^'^'2 
for small k. This is a fortunate situation, since it can be 
easily resolved by a careful calculation. To this aim, an 
explicit analytical expression (|10ap with its line-by-linc 



derivation was made publicly available [3a^. The sym- 
metry arguments of Ref. [l^ were analyzed in Ref. [2^ 
where it was shown that the symmetry argument appeals 
to a local reference system in which z-axis follows the vor- 
tex line curved by long Kelvin waves, while the Hamilto- 
nian and the interaction amplitude must be written in the 
global reference system with the z-axis oriented along the 
straight (unperturbed) vortex line. In the other words, 
the symmetry arguments are irrelevant to the problem 
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under consideration [26|, and see also [33| and 

We thus feel confident to proceed further in studying 
Kelvin wave turbulence in the framework of the Hamilto- 
nian ([6]) with the interaction amplitude (|10ap . The con- 
troversy can be reopened only when someone will present 
a checkable derivation of an alternative expression for 



3,4,5 
fc,l,2' 



which would be different from Eq. (|10aP . 



C. Some caveats and warnings 

The reader of this paper should not conclude that the 
problem of Kelvin waves in superfluids is solved. Rather, 
this is a complicated and intriguing phenomenon and 
a lot of actual problems in this field are still open or 
hardly studied. First, we stress that Kelvin waves in 
real superfluid turbulence propagate along dynamically 
bent vortex lines in a vortex tangle. Although in the 
limit of short Kelvin waves the time dependence of the 
basic vortex configuration and their local curvature can 
be neglected, longer waves are also important in real ex- 
perimental situations. Secondly and more importantly, 
even short waves could travel long distances, whereas the 
classical wave turbulence theory assumes that the wave 
"mean free path" is much less that the system size (the 
classical-quantum crossover in our case). How important 
such a finite-size is not yet understood. Numerical simu- 
lation of the full Biot-Savart equation which could clarify 
these issues is still in its infancy [3l| . 

In this paper, our analytical theory will assume that 
these limits of short and short-correlated Kelvin waves 
arc achieved. In addition we will assume that the ampli- 
tudes of Kelvin waves arc small. Then the Hamiltonian 
formulation of the problem (|6ll0ap is justified. Neverthe- 
less it does not mean that spectrum (j9bp is the only pos- 
sible solution for Kelvin wave turbulence. It known that 
approximations of weak wave turbulence leading from the 
dynamical Hamiltonian equations to the kinetic equations 
in one dimensional media are very delicate. Therefore 
in the second part of this paper we will present direct 
numerical simulations of Kelvin wave turbulence in the 
Hamiltonian formulation (|6ll0ap . This simulation allows 
us to clarify another important issue: how long should 
the vortex line be in comparison to the Kelvin wave- 
length to achieve the limit of unlimited (in the physical 
space) medium used in the theoretical analysis below. 
Also, it allows to study the case when this limit is not 
achieved and the finite-size effects are important - this is 
the case of so-called mesoscopic wave turbulence. 
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In Rcf. [l| it was assumed that the kinetic equation in 
its continuous-media hmit is apphcable. The analysis of 
its cohision integral [H, [l^l revealed that on the spectrum 
(|9ap the leading contribution appears from those terms in 
which two wave- vectors are much smaller that the other 
four and a conjecture was made (not a proof!) that the 
same kind of nonlocality should hold generally for most 
other Kelvin wave spectra, in both steady and evolving 
turbulent states. This allowed for the development of an 
effective collision integral dominated by four wave mix- 
ing of the (1 -f^ 3) type resulting in a new four-wave 
kinetic equation. The remaining two small wave- vectors 
describe the chaotic bending of the vortex lines with a 
characteristic curvature radius of the order of the inter- 
vortex distance. We should note however that even in 
the framework of the effective four-wave kinetic equa- 
tions [ij some questions have remained open. Indeed, 
the spectrum (j9b| was proven to be a valid solution of 
this equation but its uniqueness in the class of scale in- 
variant solutions and its stability have remained unclear 
and the pre-factor C^^ has not been found. In Section HIl 
we will demonstrate that the scale invariant energy spec- 
trum (j9bp exists, it is unique with a regular pre-factor, 
C^^ ~ 0.304. Moreover, our numerics of the dynamical 
equation show that the (j9bp is indeed an attracting so- 
lution for the Kelvin wave systems. Not only we observe 
the correct value of the exponent for this spectrum, but 



also a relatively good agreement for the pre-factor. This 
is remarkable because this is probably the first example 
in wave turbulence where such agreement for the pre- 
factor value has been reported in numerical simulations. 
Moreover, these numerical results demonstrate that the 
type of nonlocality conjectured for the Kelvin waves (two 
small wave vectors in a typical interacting sextet) is ro- 
bust and maintained in the attracting steady state. 

In the next section we will prove that the non-local 
theory is correct and is the one that should be used in 
future theories of quantum turbulence. After this we will 
present numerical simulations that are based on the in- 
teraction term (jlOap and which not only confirm our an- 
alytical calculations but also validate the applicability of 
the kinetic equation. These numerics also allow to study 
the regimes which are hard to treat analytically when 
the kinetic equation description breaks down due to the 
finite-size effects - so-called mesoscopic wave turbulence. 



II. ENERGY SPECTRUM OF KELVIN- WAVE 
WEAK TURBULENCE 

In this section we calculate an exact solution to the 
new effective kinetic equation. In [l[ it was shown that 
the collision integral takes the form 



St(fc,{nfc-(t)}) = — HlAk^dk2dk;[W^ 
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where the effective (1 -H- 3) interaction amplitude is 



1,2,3 _ 
fc 



kkik2kj, 



(12) 



The dimensionless parameter ^ was defined in Eq. (j9bp . 

One obvious steady-state solution of Eq. (j7bp is Uk (x 
(wfc)"^ which corresponds to the equilibrium Rayleigh- 
Jcans equipartition of the energy. A second, non- 
equilibrium solution, can be obtained by analyzing the 
equation of energy conservation, 

' = , (13) 



dt dk 
the energy flux is given by [4! 



Stfc' LUk'dk' . 



(14) 



k'<k 



In Ref. [H it was found that the kinetic equation has a 
scale-invariant solution 



Nk = Ak- 



(15) 



The scaling exponent x = 11/3 was found by power 
counting in the second Eq. (|14l) under the assumption 
that all the integral in the collision term converge. Sub- 
stituting this solution into the collision integral it was 
demonstrated in Ref. [l[ that the integral indeed con- 
verges. Finally the suggested solution has the form [l[ 
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where the dimensionless coefficient C^fj remained unde- 
termined. 

Let us now demonstrate analytically that Eq. (jl6p is 
indeed a solution of the kinetic equation. Moreover, we 
demonstrate numerically that this solution is the unique 
scale-invariant solution with a constant energy flux, and 
finally we compute the numerical value of the universal 
constant (7^^ • The first task can be achieved by con- 
sidering the second term in Eq. (|lip as a sum of three 
equivalent terms. In the first we make the Zakharov- 
Kraichnan conformal transformation jl], Q from ki,k2, k^ 
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to k[,k'2, fcg where 



kk2 



kkj 
ki 



(17) 



kj we find 



Relabeling then the dummy variable k'^ — > 
that this term, up to a factor (fci/fc)^ with ^ = 3x — 9, 
coincides with the first term in Eq. pT|) . In the second 
of the three terms, we perform the same transformation 
replacing fci fc2 , and in the third term we do the same 
but replacing ki o k^. Therefore, the entire collision 
term in Eq. (jlip can be represented as the first term 
multiplied by 



Note that the region of integration becomes the same 
for all the terms too. Due to the existence of the delta- 
function of frequencies we conclude that the integrand in 
the collision term vanishes if f = 2. This condition is 
equivalent to x = xq = 11/3. 

Next we demonstrate the uniqueness of this solu- 
tion. To do so numerically it is advantageous to non- 
dimcnsionalize the collision integral for arbitrary values 
of X as follows. 



Stk 



I{x) k 



8-3£c 



(19) 



[1 - (fci/fc)« - {k2/kf - [k^/k)^ 



(18) where a = kA/Att. Using qj = kj/k we obtain 



^ III dgidq2dq3 qi - qf) {I - qf - - 93 ) (919293)'"" S {1 - q, - q2 - Qs) S {l - ql - g| - qi) 




FIG. 1: The collision integral as a function of the scaling 
exponent. We observe that it vanishes at a; = 11/3 and that 
this cancellation is unique within the window of locality. The 
insert is just a blow-up of the neighborhood of x = 11/3. One 
clearly sees that the solution is unique. 



Combining this result with the dispersion relation, we 
finally obtain the energy spectrum 

= *^LN ^2/3 fc5/3 ■ (2^) 



This integral was computed numerically as a function 
of X, and plotted in Fig. [1] for 2 < x < 9/2 where this 
integral converges. We see that throughout the window 
of locality (where the integral converges) there is no other 
solution. 

Finally, to compute the universal coefficient Cs we re- 
turn to Eq. dU]) with Nk given by Eq. p5|) with an 
arbitrary value of x. Integrating with respect to k' leads 
to an explicit expression for the dissipation rate 



1287r 3x- ll' 



„ll-3a; 



(20) 



Using L'Hopital rule to deal with the indeterminate ratio 
of zero by zero we can rewrite this expression for x — 



The calculation of the energy spectrum for non-local 
Kelvin wave interactions constitutes the central result of 
this section. We point out that the exponent with 5/3 « 
1.67 is quite close to that of the spectrum ([5a|) with 7/5 « 
1.4. This has made it difficult for numerical simulations 
to clearly resolve the difference between both theories 
and as a result has contributed to adding more fuel to 
the controversy [l^ [2lj. In fact, a recent numerical 
simulation using a scale separation scheme for the Biot- 
Savart equation produced results consistent with the 
spectrum (Pa|) [30| . However, we believe that the scale 
separation scheme - which expands the Biot-Savart equa- 
tion in terms of non-local contributions - artificially 
damps the essential non-local behavior of the Biot-Savart 
equation, thus inducing local Kelvin wave interaction and 
the KS spectrum. 
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III. DIRECT NUMERICAL SIMULATION OF 
THE DYNAMICAL EQUATION FOR THE 
KELVIN WAVES 

A. Dynamical model: Local Nonlinear Equation 

In this section, we present a new simple model that 
contains all the necessary physics of the Biot-Savart equa- 
tion without its numerical complexity. As mentioned 
above, an exact analytical expansion of the Hamilto- 
nian ^ was carried out in [T^l leading the authors to 
the interaction term presented in Eq. (|10ap . This term is 
the leading order asymptotic of the Biot-Savart interac- 
tion in the limit when at least one Kelvin wave is much 
longer than the shortest wave in each interacting sextet. 
Obviously, the dynamically important nonlocal sextets 
in which two waves are long are also described by this 
asymptotic expression. When expressed in the physical 
space, the interaction term (jlOap leads to the Local Non- 
linear equation (LNE) which is given by 



.dw 



dt 



K d 
Att dz 



[(-; 


dw 


'] 


dw 


dz 




9^ 



(24) 



The LNE can be represented as a Hamiltonian system 
of the form ^ , with Hamiltonian which in terms of the 
physical space amplitude w{z,t) ~ x{z,t) + iy(z,t) is: 



n 



-1 


(■ 


dw 


2 


1 


dw 




AttJ 


dz 
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dz 





Nonlinear dissipation and applicability of weak 
turbulence 



To appreciate the conditions under which we could 
expect realization of the weak wave turbulence regime 
in our numerics, we should estimate the degree of the 
nonlinear dissipation given by the nonlinear resonance 
broadening parameter r^. Indeed, the kinetic equation 
for weak wave turbulence could only be valid when the 
nonlinear resonance broadening Tk is much less than the 
linear frequency Wfc. On the other hand, Tk itself could 
be estimated from the kinetic equation's collision term 
(jlip . namely from its part proportional to n^. Impor- 
tantly, even though the full collision integral converges 
on the LN spectrum, its separate parts diverge In 
particular, the part corresponding to Tk will scale as (see 
Eqs.(19) of Ref. Q)) TkUk oc rifefc giving 

Tfe (X fc . (26) 

Thus for the nonlinearity parameter we get 

— cx fc-i . (27) 



Eq. ([27)) states that the nonlinearity parameter decreases 
toward large fc, i.e. that wave turbulence becomes weaker 
as one propagates along the energy cascade to high wave 
numbers. 

On the other hand, for validity of the the kinetic equa- 
tion it is also necessary that the nonlinear resonance 
broadening is much greater than the grid spacing be- 
tween the wave frequencies Aw, i.e. Aw <C r^. If the 
wave turbulence becomes too weak, we are in danger of 
loosing wave resonances due to the sparsity (discreteness) 
of the allowed wave numbers and frequencies in bounded 
domains. If this happens, we expect regimes of discrete 
or mesoscopic wave turbulence [a Issl. [371 - 14 1| . Mesoscopic 
wave turbulence occurs when the nonlinear resonance 
broadening Tk becomes of the order of the grid spacing 
between frequencies Aw, i.e. ~ Aw. In such a situa- 
tion, discreteness effects become apparent and the wave 
amplitudes become too weak to sustain a continuous cas- 
cade. Instead, we observe sandpile behaviour [ssj . where 
energy accumulates at specific scales until it is reaches 
sufficient amplitudes to cascade. This process is then re- 
peated at smaller and smaller scales, until energy eventu- 
ally reaches the dissipation scale. When Tk <C Aw, then 
we are in a discrete wave turbulence regime, where only 
exact wave resonances occur. Resonant waves form clus- 
ters, some of which disjoint from the rest of the modes, 
and the cascade scenario is even more inhibited. 

The frequency grid spacing can be determined via 
Aw = {duJk/dk)Ak = (AKfc/27r)Afc, where Ak = 2tt/£ 
is the wave number grid spacing. Taking into account 
one can see that the ratio of Ffc/Aw is predicted to 
be independent of k implying that the wave turbulence 
description should break down when the wave modes are 
too weak simultaneously and uniformly in the whole of 
the incrtial range of k. 



C. Numerical setup 

We perform numerical simulations of the LNE to verify 
the power law scaling and coefficient of the energy spec- 
trum for Kelvin wave turbulence in a statistically steady 
state corresponding to an energy cascade from long to 
short wave scales. To achieve this, we numerically solve 
the non-dimensionalized version of Eq. p4)) in the pres- 
ence of forcing and dissipation: 



dw 



d_ 

dz 



[(-; 


dw 




dw 


dz 


'] 


'di 



iFk - iDk , (28) 



with additive forcing, -l-i-Ffc, and dissipation —iDk- The 
non-dimensionalization we use is 



Wk{t) cxp(ikz 



A 1 /4 - 47r ~ _ , 

w ~ A ' w, t ^ — -t, z = z, w 
kA 

(29) 

We implement a pseudo-spectral method with 2^^ spa- 
tial points, where the nonlinear term is fully de-aliased 
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Sim A 


Sim B 


Sim C 




610 


488 


0.6103 


A 


3.2 


2.0 


0.23 


?*: 


4.6 


4.8 


0.55 




0.35 


0.25 


0.13 


\'Hc>\/'H2 


0.25 


0.20 


0.003 -0.15 




0.3477 


0.2024 


0.031 



10' 



TABLE I: We present the non-dimensional values of the am- 
phtude of the wave action spectrum B, energy flux ik, the 
ratio of the nonlinear and linear energies \'He\/'H2 and the 



with respect to its quintic nonlincarity. We utilize 
a fourth order Runge-Kutta time stepping method to 
progress in time. The forcing profile we use is defined 
in wave number space, as is given by 



Af exp{i9k) 




if 8 < A: < 16 
otherwise , 



(30) 



where Af is a constant forcing amplitude and 0k is a 
random variable sampled from a uniform distribution on 
[0, 2n) for each k and at each time step. We add dissi- 
pation at both ends of the wave number space to enable 
the formation of a statistically steady state and to avoid 
any bottleneck effects. We use a hyper-viscosity term 
that dissipates mainly at high wave numbers and a low 
wavenumber friction on the first six wave modes. The 
dissipation profile is defined as 



Dk = I^fric H{6 ~ \k\) Wk + I^hypc 



(31) 



where, z^fric = 2.0 and fhypor = 2 x 10"^ are dissipation 
coefficients and H{-) is a Heaviside function. 

We perform three simulations, A, B and C, which 
have different forcing amplitudes (with the rest of the 
parameters being identical). The forcing amplitudes of 
the simulations are given in Table HI We allow the sim- 
ulation to evolve until we have reached a statistically 
steady state, judged by observing stationarity of the to- 
tal energy, whereupon we average over a time window to 
achieve the desired statistics. 



D. Numerical results 

For simulation A, we show in Fig. [5] the compensated 
(by k^-^/^) wave action spectrum, and compare against 
both the predictions of Eq. ([5a|) (or, in terms of the wave 
action, Nk ~ fc~^^/^) and Eq. ([^a)) (corresponding to 
Nk ~ fc^^^/^, see (US]))- The spectrum has been aver- 
aged over a long time window once a statistically steady 
state is reached. We observe a remarkable agreement 
with spectrum suggested by Lvov and Nazarenko, ([TS]) . 
for about a decade in wave number space. Moreover, 



10" 



1^ 

II 



10" 




FIG. 2: Averaged wave action spectrum Nk in simulation A, 
compensated by k^^^^ . We overlay the theoretical predictions 
of the Kozik-Svistunov and Lvov-Nazarenko spectra. 



we clearly observe a deviation from the slope of Kozik- 
Svistunov spectrum, Nk fc~^^/^. 

Also for simulation A, in Fig [3] we plot the averaged 
non-dimensionalized fiux 62 (fc) of the quadratic part of 
the energy, 'H2 ■ We observe a constant flux over a region 
of fc-space to the right of the forcing scale in agreement 
with the constant fiux scale-invariant solutions of wave 
turbulence theory. However, this flux falls off to a nega- 
tive constant at high wave numbers, which we believe is 
down to the fact that numerically we can only measure 
the ffux of the quadratic energy, which represents the 
dominant part of the total energy only for weakly non- 
linear waves. Thus, the deviation of e2(fc) from constant 
at large k can be interpreted as an increase of nonlincar- 
ity at these scales. Note, however, that theoretically one 
predicts a decrease of nonlinearity at high k, see Eq. (P?]) . 
which seems at odds with our numerical results - some- 
thing that remains to be understood. 

To add to this puzzle, we present a plot of the combined 
Fourier transform in both coordinate x and time t on the 
(fc,cj)-plane in Figured In addition to weakly nonlinear 
waves, which are narrowly distributed around the disper- 
sion relation cj = fc^, we observe a wide distribution cen- 
tered at the zero frequency. The relative intensity of such 
a "condensate" is higher at large k's: its peak is about 
1/3 of the wave intensity at the lower end of the inertial 
interval (fc = 50) and it becomes as strong as the wave 
intensity peak deeper in the inertial range (fc = 100). 
Again, we can offer no explanation for such a conden- 
sate (in the direct cascade range!) and to the fact why, 
in spite of the strong condensate, we perfectly observe 
the theoretical scaling of the Lvov-Nazarenko spectrum 
in the same range of wave numbers. 

From the numerical data of Figs. [2] and [31 and with the 
addition of the numerical measurement of ^ in Eq. (|9b[) 
we can compute the numerical prefactor C^^ of spec- 
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FIG. 3: Averaged non-dimensionalized flux of the quadratic 
energy e2{k) in simulation A. We observe a region where the 
energy flux is a positive constant, right of the forcing scale. 
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FIG. 4: The (fc, aj)-plot of simulation A representing a nor- 
malized (separately for each k) double Fourier transform of 
the wave amplitude 'w{z,t). The theoretical linear dispersion 
relation, lj — , is superimposed by the green dashed line. 



trum (|16p . where the formula for C^^ expressed in terms 
of the non-dimensional parametrization is 

C,. = {i7T)-'/'l-'/'^'^'A . (32) 

The results from the simulations and the value of the 
numerical estimate of constant C^^ are given in Table HI 

We observe remarkable agreement in the value of C^^ 
in Eq. [1^] to the numerical data from simulation B; sec 
Table [J Note that this is probably the first ever ex- 
ample of a wave turbulent system where the spectrum 
pre-factor obtained in numerics coincides with the one of 
the theoretical prediction. Indeed, realizing a pure wave 
turbulence state in numerical simulations is notoriously 
difficult because, on one hand, the waves must be weakly 
nonlinear and, on the other hand, they should not be too 
weak to overcome the finite size effects; see Sec. IIIIBI 



These effects become noticeable in our simulation A and 
they become very clear in the simulation C correspond- 
ing to the lowest forcing out of the three runs. 

In simulation B, with the forcing somewhat lower than 
in simulation A, the spectrum and the flux are very sim- 
ilar to the ones shown in Figs. [2] and [3] The spectrum ex- 
hibits a scaling range with —11/3 exponent and is clearly 
different from the —17/5 prediction of Kozik and Svis- 
tunov. However, the value of the constant C^^ is slightly 
less than the theoretical one; see TableH] This is probably 
due to the finite-size effects which manifest themselves in 
losses of wave resonances due to the discreteness of the 
fc-space when the nonlinear frequency broadening be- 
comes of the same order or less than the fc-grid spacing. 

Such an effect is most prominent in our weak-forcing 
simulation C where the observed constant C^^^ is an or- 
der of magnitude smaller; see Table ID In this run, we 
observed 'sandpile' or 'bursty' behaviour characterized 
by sudden spikes in the energies shown in Fig. O We 
observe that all energies have become statistically sta- 
tionary on the whole, however there is a clear bistable 
behaviour in the random changes to the energy that are 
intermittent in time. In fact, the ratio of the two energies 
fluctuate between \'Hq\I'H2 ^ 0.003 - 0.15. 

Obviously, the behaviour in simulation C is beyond 
the conditions which could be described by the standard 
kinetic equation. A qualitative description of this regime 
was presented in [1, IH, [s^] where the observed sandpile 
behaviour was originally conjectured. However, presently 
there is no rigorous theory describing this regime. 



E. Mesoscopic wave turbulence 

The sandpile behavior of simulation C can also be seen 
in the evolution of the wave action spectrum Nk pre- 
sented in Fig. [6] which shows two snapshots of the wave 
action spectrum at two times corresponding to a peak 
(t = 84) and trough (t — 100) of the energy. In addition, 
in Fig. [S]we include the fully averaged spectrum over the 
window 84 < T < 100 for comparison. In simulation C, 
we observe that the wave action spectrum Nk oscillates 
between the two extreme states shown in Fig. [5] at the 
same periods in which the energy fluctuations of peaks 
and troughs are seen in Fig. [51 The spectrum fluctuates 
predominately in the high wave number region. 

In Fig. [71 we present a numerical (fe,a;)-plot for sim- 
ulation C, which is similar to the one discussed before 
for simulation A, Fig. [H We sec a qualitatively similar 
picture as before, but with a more narrow wave distribu- 
tion and a weaker condensate around ui = line, which 
is natural considering the lower level of nonlinearity in 
simulation C. 

Using Fig. [3 we can find the frequency resonance 
broadening F^. This gives us a direct check if we 
are in a kinetic or mesoscopic wave turbulence regime 
d, m, [13, HI]. Numerically, we compute F^ by calcu- 
lating the width of the interval along uj in which the 
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FIG. 5: Magnitudes of the linear, 9^2, nonlinear, "He and total, 
T-i energies of simulation C at late times in units of the linear 
timescale of the forcing mode: r = Stt'^/kA/c^ where fc/ f» 10. 
We observe a statistically non-equilibrium stationary state of 
the energy. However, we also observe sharp jumps, peaks and 
troughs in the energy values. 
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FIG. 6: Two wave action spectra at times r = 84 and r = 
100, corresponding to a local maximum (red solid line) and a 
local minimum (green dashed line) in the energy respectively 
from simulation C. Moreover, we compare these wave action 
spectra with that of the fully averaged wave action spectrum 
of Fig. [3 (blue dotted line). 



dispersion relation in Fig. [7] is within half of maximum 
wave intensity at each wave number fc. We plot the result 
in Fig. [51 along with the kinetic wave turbulence thresh- 
old Ao; and the limit of the resolution of the Fourier 
transform (= 27r/r, where T is the time window of the 
Fourier transform). We observe in Fig. |5] that the fre- 
quency resonance broadening in the numerical simula- 
tion is larger than the resolution threshold, i.e. that the 
presented plot does indeed provide a useful information 
about Ffc. Further, we see that F^. grows to higher fc con- 
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FIG. 7: The (fe,a;)-plot of simulation C. The random waves 
are characterized by the black curve, which is a normalized 
double Fourier transform of the wave amplitude w{z,t). The 
theoretical linear dispersion relation, oJk = k^, is superim- 
posed by the green dashed line. 




FIG. 8: We plot the nonlinear frequency broadening Tk of 
the (fe, tj)-plot from simulation C in the red solid line. We 
overlay the resolution of the fast Fourier transform (green 
dashed line) and width of the frequency spacing Auj (blue 
dotted line). The black dashed line corresponds to a linear fit 
to Ffc. 



sistently with theoretically predicted linear dependence 
(PS]) . Finally, we see that the values of F^ are about fac- 
tor three above the line corresponding to the threshold 
of mesoscopic wave turbulence F;; ~ Alo. Since numer- 
ical factor of order unity are not controlled by such a 
simple estimate, we can conclude that such values of F^ 
are consistent with the condition of the mesoscopic wave 
turbulence regime. 

Finally, it is also meaningful to ask if any mesoscopic 
effects, e.g. the sandpile behaviour, can be seen in the 
physical space. We observe an interesting two-scale pic- 
ture, which is consistent with the scale separation as- 
sumed in the non-local theory. Further, we observe that 
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FIG. 9: We plot the physical space intensity I{z, t) — 
1^(2,^)1^ for simulation C during a time period where the 
energy is fluctuating. Here black implies strong intensity and 
white implies low intensity. The large scale structures ob- 
served correspond to the forcing scale. 



intensity of the slow (nearly horizontal) large-scale rolls 
suddenly drops at the moments corresponding to the 
drop of the total energy in Fig. [5] marking onsets of 
the "sandpile tip-overs". The same large-scale rolls are 
the main contributors into the low-frequency condensate 
component. 



IV. CONCLUSION 

In conclusion, we have shown that the spectrum pro- 
posed in [l| for Kelvin wave turbulence is the exact and 
unique solution of the kinetic equation and therefore is 



the one that should be physically realized. All of the 
theoretical predictions have been confirmed via numer- 
ical simulations of the LNE, including both the scaling 
exponent and the spectrum pre-factor. We stress that 
the LNE is only applicable for nonlocal Kelvin wave tur- 
bulence. Thus, our numerical result that the (theoreti- 
cally predicted) nonlocal spectrum is an attracting state 
demonstrates that the nonlocal dynamics is self-sustained 
and robust and the LNE model is self-consistent. 

We also studied a mesoscopic regime of Kelvin wave 
turbulence which occurs at low amplitudes. We discov- 
ered that the system suffers from intermittent energy 
bursts that produce fluctuations in the wave action spec- 
trum at high wave numbers. This agrees with the sand- 
pile scenario of evolution which was predicted for meso- 
scopic wave turbulence in d, [H, 113, HI] . 

Finally, numerically we observed some effects which 
we cannot presently explain. Namely, in addition to 
the weak dispersive waves we have found a strong low- 
frequency condensate component over a large range of 
scales covering the direct cascade range. In the small- 
scale part of the inertial range the condensate is of a 
similar strength to the wave component, and yet, some- 
how, it does not alter the theoretical (non-local) weak 
turbulence scaling. 
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